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We report the statistical properties of spherical steel particles rolling on an inclined surface being 
driven by an oscillating wall. Strong dissipation occurs due to collisions between the particles and 
rolling and can be tuned by changing the number density. The velocities of the particles are observed 
to be correlated over large distances comparable to the system size. The distribution of velocities 
deviates strongly from a Gaussian. The degree of the deviation, as measured by the kurtosis of 
the distribution, is observed to be as much as four times the value corresponding to a Gaussian, 
signaling a significant breakdown of the assumption of negligible velocity correlations in a granular 
system. 

PACS number(s): 81.05.Rm, 05.20.Dd, 45.05,+x 
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From the study of landslides to the understanding of 
the formation of galaxies, systems of macroscopic parti- 
cles are of fundamental interest for physicists and engi- 
neers . To treat a granular system statistically, tem- 
perature must be replaced by a more suitable quantity. 
Granular temperature has gained considerable attention 
over the past two decades as a useful quantity for describ- 
ing the properties of rapid granular flow (^-j^] . If the colli- 
sions between particles are elastic, then we would expect 
on very general grounds that the velocity distribution is 
Gaussian and the velocity components are uncorrelated. 
However, dissipation due to inelasticity and friction can 
introduce correlations that makes these properties of the 
velocity distribution suspect [ pp[ . For example, a simu- 
lation of a system of inelastic particles yields long-range 
velocity correlations 0]. 

The first experiments with steel particles excited by a 
vertically vibrated container indicated that the velocity 
distribution is Gaussian || . Using advances in high speed 
digital imaging, later experiments have shown deviations 
in the tails of the distribution from a Gaussian [p|-pr|. 
It also has been claimed that the distribution of the 
velocity component, v x , in the direction perpendicular 
to the direction of the vibration, may be described by 
the functional form R(v x ) = exp[— (\v x \/ a x ) 3 / 2 ], where 

a x = (w^) 1 ^ 2 [ p]JT^1 . However, these experiments were 
conducted with particles that are relatively elastic, and 
the velocity correlations were either not measured or 
claimed to be negligible. Particles that collide while 
rolling lose more energy in comparison to particles that 
collide in free space, due to frustration at impact. In ear- 
lier work |H| , this effect lead to strong deviations even in 
a dilute system where the mean free path was of the order 
of the dimension of the container. However, dissipation 
at the boundary was significant and the distribution did 
not represent the bulk properties of the particles. These 
observations lead us to ask if there are correlations in 
dense systems where the mean free path is much smaller 



than the system size, and if these correlations affect the 
distribution of the velocities. 

In this Letter, we address these questions by measur- 
ing the velocities of particles rolling on an inclined plane. 
This geometry slows the dynamics and allows us to vary 
the dissipation over a wide range. Energy is supplied to 
the system by means of an oscillating wall. We report 
the density and velocity distributions as a function of 
the number of particles N in a regime where the distri- 
butions are insensitive to the frequency or phase of the 
drive. The dissipation increases with N because the time 
between collisions as well as the mean free path decreases. 
Because the density of the particles becomes higher near 
the driving wall as N is increased, our analysis of the 
velocities is done in a narrow region where the density 
and the velocity fluctuations are relatively uniform. 

We observe that significant correlations in the veloci- 
ties of the particles that increase in strength and range 
as N increases are found. The distribution, P(v x ), of 
the velocity components of the particles perpendicular 
to the direction of the oscillation deviates strongly from 
a Gaussian. At low N , corresponding to lower dissipa- 
tion, P{v x ) may be fitted to R{v x ) as in Refs. |Il] , |l2| ], but 
P(v x ) deviates strongly from this form as N is increased. 
The kurtosis of P(v x ) is observed to be as much as four 
times the value corresponding to a Gaussian distribution, 
signaling a significant breakdown of the assumptions of 
negligible velocity correlations in our system. 

The experimental system consists of stainless steel 
spheres with diameter d = 0.3175 cm rolling on a 
29.5 cm x 27.0 cm glass plane that can be tilted from 
= < 8°; a schematic of the experimental setup is 
given in Ref. [fL3|| . The coefficient of rolling friction is 
observed to be less that 1 x 10~ 3 , and therefore energy 
loss during rolling is negligible compared to the loss that 
occurs due to a collision |l4|,[l5| • The number of particles 
is varied from N — 100 — 1000. To track the positions of 
the particles we utilize a high speed digital camera (Ko- 
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dak Motioncorder) at a frame rate of 250 per second. A 
centroid technique is used to follow the particles at sub- 
pixel accuracy, allowing us to resolve positions to within 
0.13mm and velocities to within 5.2 x 10 _4 cms — ^ The 
statistical properties are obtained by averaging over time 
using at least 8500 frames. 
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FIG. 1. (a) The average number density, p(y), in the 
{/-direction for various values of N; 1000 (•), 750 (o), 500 (o), 
200 (A), 100 (+). The driving wall is at y = 26 cm. (b) The 
radial distribution function g(r) for the same values of N at 
/ = 6 Hz (d is the particle diameter). The plots of g(r) have 
been shifted for clarity, and the dashed lines show g(r) = 1. 
To show that the details of the driving do not influence the 
correlation of the particles, g(r) for N = 750 at / = 6 Hz (o) 
and / = 20 Hz (o) are shown. 

Energy is added to the system by means of the bottom 
wall which is attached to a solenoid powered by a wave- 
form generator. The driving signal is not sinusoidal, but 
a periodic pulse with an amplitude A ~ 2.5 d and a veloc- 
ity Uwaii ~ 20 cm s -1 . When the system is tilted past 9 = 
2.0° and the driving frequency / > 4.0 Hz, the positions 
of the particles are observed to be independent of the 
phase of the drive. Below 4.0 Hz and for N > 200, we ob- 
serve particles forming subharmonic patterns similar to 
those observed in vertical vibrated quasi-two-dimensional 
geometries |M . The nature of the velocity fields and den- 
sity distribution in this latter regime will be discussed 



elsewhere [[L7| . All of our results reported here are at a 
fixed 6 = 3.0°, A = 2.5 d, and / = 6.0Hz. For these pa- 
rameters phase dependence is absent in the density and 
velocity distributions. 




FIG. 2. The spatial correlation function of the velocities, 
C||(r), as a function of separation r/d, for N = 1000 (•), 750 
(o), 500 (o), and 200 (A). Inset: To show the nonexponen- 
tial decay of Cu(r), we plot In Cm (r) versus In (r/d). The solid 
lines are least squares fits to the data. The slope a = —1.23 
(•), -1.38 (o), -1.55 (o), -1.90(A). The statistical error bars 
are smaller then the symbols. 

The x and y axes are chosen to be along the driving 
wall and perpendicular to the driving wall respectively. 
The origin corresponds to the point where the top wall 
intersects the side wall. The distribution of the particles 
in the x-direction is uniform. The y-dependence of the 
number density, p(y), is obtained by averaging the num- 
ber of particles in bins of width d and dividing by the 
area of each bin (see Fig. |l|(a)). The density p(y) is not 
constant because of gravity and the action of the driving 
wall, although for systems with N > 100, p(y) shows a 
plateau where there are very small density gradients, sim- 
ilar to earlier observations Jl8],|l9| . Most of the particles 
remain suspended above the driving wall and receive ran- 
dom kicks from the particles that are between the bulk 
and the driving wall. 

To ensure that the deviations in density over the com- 
pacted region do not influence the distribution of veloc- 
ities, we varied the region over which the analysis of the 
velocities was performed. We found that the distribu- 
tion of velocities is independent of the distance y in the 
compacted region. Therefore, we measured the velocity 
distributions in a region of width bd beginning 8c? away 
from the driving wall. 

To further characterize the compacted region in our 
experiment, we measure the radial distribution function 
g(r) as a function of the distance r between the centers 
of the particles M . Our results for various N are shown 
in Fig. ^(b). The plots demonstrate a transition from a 
gas- like state for TV < 200, to a more liquid-like state for 
N > 500. 
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We now discuss the detailed properties of the velocities 
of the particles. The spatial correlation function of the 
velocities is given by 



where i and j label particles separated by a distance r. 
With this definition, two particles with parallel (antipar- 
allel) velocities give a correlation of +1 (—1). We find 
that long-range velocity correlations are present for all 
N as shown in Fig. |[ We have plotted lnCii(r) versus 
In (r/d) in the inset of Fig. ||. The r-dependence is non- 
exponential, and least squares fits to In C\\(r) suggest a 
possible power law decay. However, our range of r is lim- 
ited by the finite size of our system. The values of the 
slope a are given in the caption of Fig. |[ 

We now discuss how these correlations affect the dis- 
tribution of velocities in the x and y directions (see Fig- 
ure |^). The velocity components v x and v y have been 

scaled by a x = (y 2 ) 1 ^ 2 and a y = {v 2 ) 1 ^ 2 ■ The maxi- 
mum of the distribution P(v x ) is scaled to be unity for 
clarity and is given in Fig. ||(a) for various N. A Gaus- 
sian fit is shown for comparison. For low densities, we 
find that the form of P(v x ) is similar to that observed 
in previous experiments with low dissipation ^,0. For 
N = 100 (+), the form of P(v x ) can be fitted by R(v x ) = 
exp[—(\v x \/<T x ) 3 / 2 ], shown by the dashed line. However, 
as N is increased, the deviations are observed to become 
stronger and cannot be described by this power law form. 

The asymmetry of P(v y ) seen in Fig. |^(b) occurs be- 
cause the particles moving toward the piston have lost 
energy because of collisions. The distributions become 
more asymmetric as N is increased because of the greater 
number of collisions between particles. To characterize 
P(v y ), we calculate its skewness (third moment), and 
find that it ranges from 1.62 - 4.45 for TV = 100 - 1000 
respectively. The effects of this asymmetry on P(v x ) can 
be determined by discriminating particles with v y < 0. 
After performing this "filtering," we find that P(v x ) is 
insensitive to the sign of v y . Thus the asymmetry of 
P(v y ) is not the cause of the deviations in P{v x ) from a 
Gaussian. 

By using the kurtosis or flatness of the distribution 
P(v x ) given by F x — (v x ) / (v 2 ) , we can quantify the 
deviation of P(v x ) from a Gaussian. If the distributions 
were Gaussian, F x = 3.0. If the distributions were given 
by R(v x ), F x = 3.762. Figure |(c) shows F x as a func- 
tion of N. The increase in F x shown in Fig. ^(a) implies 
that P(v x ) deviates more strongly from a Gaussian as 
dissipation is increased. The linear dependence of F x on 
N may be due to the finite range of N. 
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FIG. 3. The distribution of velocities for N = 1000 (•), 
750 (o), 500 (o), 200 (A), and 100 (+). (a) The veloc- 
ity distribution P(v x ) for the region described in the text. 
Each distribution is normalized to be unity at its maximum. 
The deviation from a Gaussian (solid line) are significant 
for all N and show a systematic increase as a function of 
N. In the most dilute case (+) P(v x ) can be described by 
R{v x ) = exp[— (\v x \/a x ) 3 ' 12 ] shown by the dashed line, (b) 
The distributions P(v y ) are observed to be asymmetric due to 
the driving wall and dissipation. The asymmetry is observed 
to increase as -/V is increased. However, this asymmetry does 
not cause the deviations from a Gaussian in P(v x ) (see text). 

These results show that there are strong correlations, 
both in the positions and in the velocities, as a conse- 
quence of the high dissipation that occurs due to colli- 
sions. If steel particles collide head-on in free space, they 
lose energy as determined by the coefficient of normal 
restitution given by r\ rj 0.93. However, when rolling 
particles collide, the translational velocity changes at the 
instant of collision, but the angular velocity does not. 
Therefore, the rolling condition is no longer satisfied and 
the particles slide for a short duration. Because the slid- 
ing friction is about 100 times larger than the rolling fric- 
tion for steel on glass, the rolling condition is established 
within a few particle diameters |T3]-|T5|]. Hence, particles 
lose more energy than they would by considering only 
the normal coefficient of restitution. Measurements for 
our system show that the effective coefficient of restitu- 
tion T]eS ~ 0.5 |L3]. Because the effect of inelasticity is 
to reduce relative motion plj, this smaller value implies 
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that the particles tend to cluster and stream when driven, 
leading to the observed spatial and velocity correlations. 

To estimate the increase in dissipation with N, we cal- 
culate the mean free path i from the relation given in 
Rcf. p2|. The results for I as a function of N are shown 
in Fig. ^(b). We find that I decreases much faster than 
the estimate based on an uniform spatial distribution of 
particles. By dividing i by the root mean square veloc- 
ity, we obtain the mean collision time between particles. 
This time is observed to decrease from 0.83 s to 0.059 s 
as N is increased. The increase in the range of the ve- 
locity correlations and the deviations from a Gaussian in 
P(v x ) is due to the higher dissipation in the system for 
larger N. As TV is increased, the collision rate increases, 
which gives rise to higher dissipation. As the latter is 
increased, the overall particle velocities decrease. This 
decrease in turn leads to more compaction and smaller 
mean free paths. 
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FIG. 4. (a) The kurtosis F x of P(v x ), as a function of N 
(F = 3 for a Gaussian.) The solid line is a least squares fit to 
the data, (b) The mean free path I as function of the number 
of particles N. (c) The standard deviations a x (o) and a y 
(a) (measure of granular temperature) as a function of N, 
taken from the distributions in Fig. ^. The difference of a x 
and o v reiterates that the granular temperature is a nonscalar 
quantity. The line is a guide for the eye. 

To ensure that the means of energy input is not effect- 
ing the system, we have varied the driving parameters. 
We find that the statistical properties do not depend 
strongly on how the system is driven above a thresh- 
old / > 4 Hz. Increasing / and the amplitude A of the 
oscillating wall by a factor of five has no effect on the 
spatial and velocity distributions for fixed N . For ex- 
ample, Fig. p](b) shows that for N = 750, g(r) remains 



unchanged if / is increased by a factor of 3.3 and A is 
increased by a factor of 3.0. 

In Fig. |4|(c) we show a x and a y as a function of N 
to confirm the anisotropy of the granular temperature 
for our system. The plot also shows that the average 
energy per particle decreases as N is increased. This de- 
crease occurs due to the higher collision rate between the 
particles that removes energy from the system through 
dissipation. 

In summary, strong velocity correlations in dense gran- 
ular gases have been experimentally reported for the first 
time. These correlations are observed to cause the veloc- 
ity distributions to deviate significantly from a Gaussian. 
Our experiments also suggest that the velocity distribu- 
tions depend on the degree of dissipation in the system 
and the form R(v x ) = exp[— (\v x \/a x ) 3 / 2 ] of the velocity 
distribution is not universal. 
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